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ABSTRACT 

Bouchet et al. (1992) showed that in an open or closed Universe with only 
pressureless matter, gravitational instability from Gaussian initial conditions 
induces a normalized skewness, S3 = (S 3 ) (S 2 ) , that has only a very weak 
dependence on the nonrelativistic-matter density. Here we generalize this result 
to a plethora of models with various contributions to the total energy density, 
including nonrelativistic matter, a cosmological constant, and other forms of 
missing energy. Our numerical results show that the skewness (and bispectrum) 
depend only very weakly (<; 2%) on the expansion history. Thus, the skewness 
and bispectrum provide a robust test of gravitational instability from Gaussian 
initial conditions, independent of the underlying cosmological model. 



Subject headings: cosmology: theory — galaxies: large scale structure of the 
universe — galaxies: clustering — galaxies: statistics 



1. INTRODUCTION 

One of the robust predictions of inflation that distinguishes it from numerous 
alternatives such as topological-defect or isocurvature models is that large-scale structure 
grew from a Gaussian distribution of primordial density perturbations. However, even 



^amionQphys. columbia.edu 
2 ari@astro. columbia.edu 



-2- 



if the primordial distribution was Gaussian, subsequent gravitational evolution would 
introduce deviations from Gaussianity. A Gaussian distribution is symmetric with 
respect to overdense and underdense regions, but as gravitational amplification of density 
perturbations proceeds, overdensities can become arbitrarily large. On the other hand, any 
given region can become only so underdense (since densities must be positive). Thus, the 
resulting distribution is necessarily skewed. 

Peebles (1980) calculated perturbatively the normalized skewness induced by 
gravitational instability of initially Gaussian perturbations in an Einstein-de Sitter 
universe. More precisely, if S(x) = [p(x) — p]/p is the fractional density perturbation 
and p(x) is the density at position x, then in an Einstein-de Sitter universe, the statistic 
S 3 = (<5 3 ) (<5 2 ) = 34/7. Bouchet et al. (1992) generalized this result to an open or 
closed universe containing only pressureless matter and found that the prediction for 53 
remains the same to < 2% for any reasonable value of Qq > 0.1. Thus, they showed that 
measurements of S3 from galaxy surveys should provide a robust test of Gaussian initial 
conditions, independent of il , but only if the universe consists of just pressureless matter. 

However, slow-roll inflation predicts that the universe is flat, while observations seem 
to point to a matter density Qo — 0.2 — 0.5 < 1. Therefore, if inflation is correct, then 
there must be some other missing energy density. The most widely considered form for this 
density is a cosmological constant (with pressure given by p = — p), but theorists have also 
realized that this missing energy might be characterized by any of a number of equations of 
state. Among these are some energy density that scales as a~ 2 (if-matter), where a(t) is 
the scale factor of the universe ( |Kolb l989| ; |Kamionkowski fc Toumbas 1996| ; Pen fc Spergel 



|1997|) , something akin to an evolving cosmological constant, possibly generated by the 
dynamics of some scalar field ( |Goble, Dodelson, fc Frieman 1997] |Silveira fc Waga 1997 



[l'urner fc White 1997| ; |Oaldwell, Dave, fc Steinhardt 199"5| ), or some combination of several 



forms of missing energy. 

To test for primordial Gaussianity of density perturbations in this cornucopia of 
missing-energy models, we calculate the induced skewness in flat models with a cosmological 
constant and/or some other form(s) of energy. For completeness we also consider open 
and closed universes with a variety of equations of state. For all cosmological models 
investigated, the numerical results for S3 never deviate from 34/7 by more than 2% for 
models with Q > 0.1, and no more than about 1% for Q > 0.3. Thus, we conclude that 
the normalized skewness does indeed provide a robust test for primordial Gaussianity, 
independent of the nature of the species that constitute the total energy density. 

The bispectrum (Fourier transform of the induced three-point correlation function) has 
also been calculated, as a byproduct of the calculation of the skewness, for Gaussian initial 
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conditions ([Peebles 1980| ; Fry 1984 ; Goroft' et al. 198"6 ) in an Einstein-de Sitter universe 



and in a universe with pressureless matter ( |Bouchet et al. 1992| ). Here we provide results 



for the bispectrum for the models we investigate. 

The details of our calculation are presented in the next Section. We write down the 
differential equation, for an arbitrary expansion history, that must be solved to obtain 
the growth of the skewness. In Section 3, we provide numerical results for the skewness 
and bispectrum for several forms of missing energy and for some models with various 
combinations of missing energy. We also provide an analytic approximation to these 
quantities as a function of Qq for flat cosmological-constant models. We then make some 
concluding remarks. 



2. CALCULATION 

The calculation of the bispectrum and skewness is outlined in detail in Peebles (1980) 
and Fry (1984), and we simply highlight the relevant steps below; Bouchet et al. (1992) 
employ an alternative approach based on the Zeldovich approximation. The Friedmann 
equation for the expansion rate as a function of redshift, z, is given by 

H(z) = ~ = H E(z); E{z) = ^ {1,(1 + ^O^d, (1) 

where H is the present-day value of the Hubble parameter, the dot denotes a derivative 
with respect to time t, and J2i ^ = 1 where fl, represents the contribution to the overall 
energy density from species i having equation of state p = Wip. In this formulation, the 
curvature of the universe contributes an amount Qk to the total energy density, and yields 
the term Qk{1 + z) 2 in the sum in equation ([!]). Thus, in a universe with, for example, 
nonrelativistic-matter density Qq and cosmological-constant contribution to closure density 
of E(z) = yQo(l + z) 3 + fl\ + (1 — Q — CI a) (I + z) 2 . The deceleration can be written 
as 

~ n = H 2 F(z), (2) 

and F(z) can be obtained by differentiating equation (p. For example, for the E(z) given 
above, F(z) = fl A - Q {1 + zf /2. 

If the fractional density perturbation is small (5 1), the equations of motion for 
<5(x, t) can be solved perturbatively using the expansion <5(x, t) = 8^(x,t) + 5^ 2 \x,t) + • • •, 
with <^ 5^K The perturbative equations for the lowest-order term, turn out to 
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be separable, 5 ( - 1 \x, t) = Si( y x)D 1 (t), and the familiar linear-theory growth factor satisfies 



a 2 V a 



£>i = 0. 



(3) 



The next-order equation of motion is 



+ 




^(A Mi )(A! 



where 



«5i(x) 



(4) 



(5) 



The homogeneous part of equation (g) is the same as that for equation (H), and it turns 
out that the solution for 5^ 2 ^ is dominated by the inhomogeneous terms on the right-hand 
side. Since the inhomogeneous part is the sum of three terms with two different time 
dependences, the solution for 6^ is not obviously separable. However, the complete solution 
is simply the sum of the solutions to the equations given by the homogeneous equation 
driven by each of these sources separately. We may thus write the solution as 



S {2) = (£> 2 ,a + D 2)b )6l + {D 2>a + 2D 2;b )6 1)i {A 1)i ) + A^A^Ai 



where the second-order growth factors satisfy 



D 2 , a + 2-D 2; , 
a 



2 



;n # 2 (-) D 2 , a 



(6) 



(7) 



D 2ib + 2-D 2fi 
a 



D 



2,6 



n 2 



with the initial conditions that the growth factors and their first time derivatives are zero 
at t — 0. Although it appears that there are two differential equations (for D 2 ^ a and D 2 ^) to 
be solved for the second-order solution, it can be verified that D 2b = (D\ — D 2>a )/2. Thus, 
there is really only one independent second-order differential equation. The first of these 
equations (|7|) can be solved analytically for the Einstein-de Sitter model with the result 
that D\ oc t 2 / 3 and D 2>a = (3/7)Df. For more general models, equations (|]) and (|7|) can be 
integrated numerically. 
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It is then straightforward to follow the arguments given in §18 of Peebles (1980) to find 
the normalized skewness. In terms of the quantity fi = D 2a jD\ [this is twice the parameter 
k in Bouchet et al. (1992)], it is given by 



S 3 = 4 + 2/i = — + - [-n 



1 



(9) 



We have thus enumerated a simple prescription for calculating S3 in a model with Gaussian 
initial conditions and an arbitrary expansion history. 

In terms of the parameter /i, the second-order correction to the density contrast is 
( IBouchet et al. 1992|) , 



5 m 



D 2 r 
^ (1 



+ fi)5 2 + 2<Ji, i Ai >i + (1 - Ai^Ai 



■'j 



(10) 



If we compare equation (11) in Fry (1984) with equation (A2) in Goroff et al. (1986), it 
is clear that the unsymmetrized scaled bispectrum for arbitrary expansion history can be 
written [cf. equation (Al) in Goroff et al. (1986)] 



P 2 W (ki,k 2 ) 



(1 + ^ 4. ki • ks fh. + ^i) + (1 _ fA ( ki • k2 



3. RESULTS 

The solid curve in Fig. |l] shows the normalized skewness ^3 as a function of the 
no nrelativistic- matter density Qq in an open, flat, or closed model containing only 
pressureless matter; i.e., in a model with E 2 (z) = f2o(i + z) 3 + (1 — fio)(l + z ) 2 - [To 
clarify, the results do not really depend on the geometry of the universe, but rather on the 
expansion history E(z).] This reproduces the result of Bouchet et al. (1992) who found 
that the quantity \i is well approximated for 0.05 < Q < 3 by 

H - ^0 V63 ■ (12) 

The short-dashed curve in Fig. [I] shows our results for the normalized skewness S3 as 
a function of the nonrelativistic-matter density Qq in a flat cosmological-constant modelQ 
that is, E 2 (z) = £l (l + zf + (1 - Q ) (note that A < for models with fi > 0). We see 



3 Bouchet et al. (1995) have also integrated these equations numerically for the cosmological-constant 
universe 
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Fig. 1. — The normalized skewness S3 as a function of the nonrelativistic- matter density 
Qo in a universe with only nonrelativistic matter (solid curve), a flat cosmological-constant 
model (short-dash curve), and a flat model with an energy density that scales as (1 + z)' 1 
(long-dash curve). 
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that the difference between S3 and 34/7 is even smaller in a cosmological-constant model 
than in an open model with the same Jl . We find that p, is well approximated in these 
models by 

p * 3 -n^. (13) 

To generalize even further, the long-dashed curve in Fig. [1| shows the normalized 
skewness S3 in a flat model with a nonrelativistic-matter density Qq plus some hypothetical 
matter density with p = — 2p/3, scaling as (l+z) -1 ; that is, E 2 {z) = Q (l+z) 3 +(l— Q )(l+z) 
(we will refer to this as "Ox-matter" below). Again, the correction to S3 is very small for 
any reasonable value of Qq. The normalized skewness is well approximated in these models 
by 

Ss - |nJ /7T . (14) 

A flat model with nonrelativistic matter and some i^-matter density, such as 
nonintersecting cosmic strings, which scales as (1 + z)~ 2 ( |Kolb I§g3 ; |Kamionkowski & 



Toumbas 1996| ; [Pen fc Spergel 199?] ), has the same E(z) as the standard open model with 



the same nonrelativistic-matter density (since wk = —1/3). Thus, the skewness will be 
given by the solid curve in Fig. [l]. 

Of course, the cases shown in the Figure do not exhaust the full range of E(z), and 
we can think of no way to systematically study all possible cases. However, we have 
tried numerous strange combinations of the types of matter considered above, and in no 
case do we find any deviation from 34/7 significantly larger than that in the open model. 
Several examples are shown in Table 1. In addition to flat Qq < 1 models with various 
combinations of cosmo logical constant, K- matter, and f^-matter (including some with a 
negative cosmological constant), we also consider a closed cosmological-constant model that 
is just consistent with quasar-lensing statistics flWhite fc Scott 1996Q and an fl = 10 closed 
model ( Harrison 1993|) . We include these models to demonstrate the robustness of S3 to 



dramatic variations in the expansion history. We do not suggest that all of these models 
are observationally tenable. These numerical experiments lead us to believe with good 
confidence that S3 must differ from 34/7 by no more than 2% in any cosmological model 
with Qq > 0.1, and to no more than ~ 1% for any model with > 0.3, if structure grew via 
gravitational instability from a Gaussian distribution of primordial density perturbations. 
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Table 1: Predicted skewnesses and fi. The second to last column lists the percentage 
deviation of S 3 from 34/7. 

4. DISCUSSION 

We have generalized the result of Bouchet et al. (1992) — that the normalized 
skewness and bispectrum are only very weakly sensitive to Qo m models with only 
pressureless matter — to general Friedmann cosmologies with arbitrary contributions to the 
energy density, including the simplest cosmological-constant model. The robustness of 
these measures allows for excellent model-independent tests of structure formation from 
gravitational instability of a Gaussian distribution of primordial perturbations. 

We have only considered the effect of the expansion on the skewness and bispectrum. 
One might wonder whether perturbations in the scalar field in models with a dynamical 
scalar field might affect the growth of structure. If so, then it is conceivable that the 
skewness would deviate from the gravitational-instability prediction. However, structure 
formation in this case would not be due simply to gravitational instability, as we have 
assumed here. Analytic approximations suggest that on small scales 100 hr x Mpc), 
scalar-field ordering should not significantly affect the skewness flJaffe 1994j) . Thus we 
expect that in (non-topological) scalar-field models, the prediction S3 ~ 34/7 should hold 
for Gaussian initial conditions. 

In practice, other factors must be considered when comparing measurements of the 
observed skewness and bispectrum with the predictions above. For example, the skewness 
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calculated above is that for an unsmoothed density field, while the observed distribution 
is intrinsically discrete. Bernardeau (1994) considered the effects of smoothing, and it 
is straightforward to apply the discussion therein to generalize our results for S3. In 
particular, for the case of smoothing with a spherical top-hat filter for scale-free power 
spectra P(k) oc k n , 

S 3>s = S 3 -(3 + n), (15) 

where 5*3 jS is the result for the smoothed distribution. In addition, the skewness and 
bispectrum presented here are those for the underlying mass distribution. Galaxy surveys 
probe the luminous-matter distribution, which if biased relative to the matter, will be 
characterized by a different skewness and bispectrum. If the galaxy fractional density 
perturbation is written as 



5 9 (x,t) = M(x,t) + |5 2 (x, £) + ■■•, (16) 



where b± is the linear bias term, 62 the first nonlinear term, and so on, then the skewness of 
the observed galaxy distribution is given by ([Fry fc Gaztanaga 1993Q 



S3,g - fr^^s + ( 17 ) 



and similar corrections have been obtained for the bispectrum ( [Fry 1994]) . More realistic 



calculations involving, for example, true scale-dependent power spectra or bias evolution 
have been obtained using numerical techniques ([Fry 1996| ; |Jing fc Borner[ |Gaztanaga 



|fc Bernardeau 1998| [Buchalter fc Kamionkowski 1998]) . The model-independence of the 



predicted skewness and bispectrum for the matter distribution makes determinations of the 
bias from measurements of these quantities ( |Matarrese, Verde, fc Heavens 1997| ; Puchalter 



fc Kamionkowski 1998| ) that much more robust. 



The robustness of the scaled skewness to the expansion history supports the notion 
that the geometry of gravitationally collapsing objects is determined almost exclusively 
by the initial conditions. This leads us to speculate that the predictions for the higher- 
order moments and correlation functions (e.g., kurtosis, four-point correlation function, 
trispectrum, and so forth) will be similarly independent of the expansion history. It would 
be interesting and straightforward to check this hypothesis numerically. 
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NAG5-3091, NSF AST94-19906, and the Alfred P. Sloan Foundation. 
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